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Abstract 

A 3 sigma excess has been recently announced by ATLAS in events 
with Z-peaked dilepton pairs, jets, and large transverse missing energy. 
We interpret this finding in the context of composite Higgs / RS theo¬ 
ries. We find that composite Higgs theories with custodial symmetry 
protection to the Zbb coupling predict a significant contribution to 
ZZbb (and to hhbb ) final states coming from heavy gluon decays to 
pairs of bottom-partner vectorlike quarks. The heavy gluon to vector¬ 
like quarks signal is largely accepted by the ATLAS selection if one 
of the Z boson in the ZZbb final state decays leptonically and the 
other to neutrinos. For a bottom partner of ~900 GeV, we find that 
the ATLAS excess can be reproduced by composite Higgs models, in 
an experimentally allowed parameter space, for heavy gluon masses 
roughly in a range 1.87 - 2.15 TeV and for heavy gluon couplings to 
light quarks within ~ (0.3—0.65)55. We briefly discuss the implication 
of this result for future experimental tests. 
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1 Introduction 


An excess above the expected number of Standard Model events, with a sta¬ 
tistical significance of three standard deviations, has been recently measured 
by ATLAS in events containing a same-flavor opposite-sign dilepton pair, 
jets, and large transverse missing energy [lj. The excess is only found in the 
on-Z region, for dilepton mass near the Z peak. 0 

The on-Z ATLAS search is focused on generalized gauge mediated (GGM) 
supersymmetry breaking models and aims at selecting a signature where pair- 
produced gluinos decay via neutralinos to a gravitino plus a Z boson. Two 
examples of GGM models are considered by ATLAS to interpret the results. 
Recent studies [d[j5] have also analyzed the ATLAS results in the context of 
supersymmetry. In |3j it is found that “only relatively light gluinos, together 
with a heavy neutralino decaying predominantly to Z boson plus a light 
gravitino could reproduce the excess”, while j4) shows that GGM models 
are unlikely to reproduce the ATLAS results and that alternative scenarios 
with two massive neutralinos, which can be possibly realized in NMSSM, 
are favored. The study in [5] arrives at a similar conclusion and states that, 
after combining the relevant constraints from LHC searches for new Physics, 
GGM models cannot explain the ATLAS excess. 

In this paper we present an interpretation of the ATLAS results in the con¬ 
text of composite Higgs 16]/ RS |7| theories. Composite Higgs models are 
compelling theories to solve the hierarchy problem and to explain the elec- 
troweak symmetry breaking. In the most natural parameter region of the 
theory |8j[9], new vector resonances are above the threshold for decays into 
new vector-like top partners and we will show that the ATLAS excess can in¬ 
deed be produced by a heavy gluon resonance decaying to vector-like quarks. 
In particular, we will show that in model with a custodial symmetry pro¬ 
tecting the Zbb coupling 10 heavy gluon decays to bottom partners can 


give a large contribution to hhbb or ZZbb final states. In the latter case, If 
one of the Z decays leptonically and the other to neutrinos, thus producing 


1 A similar analysis in the same channel considered by ATLAS has been also performed 
by CMS [2]. Probably due to the different kinematic selection, as also speculated in 1 
(in particular CMS does not apply a cut in H T ), CMS does not find any excess of events 
above the background in the on-Z region. On the other hand CMS claims a 2.6 a excess 
in the off-Z region, where no deviation from SM expectations are observed by ATLAS. In 
our study we will not try to explain the CMS excess and we will not analyze the CMS 
results. 
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large missing energy, the heavy gluon to vector-like quarks signal is largely 
accepted by the ATLAS selection [l]. 

We will analyze in detail the heavy gluon to bottom partners to ZZbb sig¬ 
nature, considering as a benchmark a concrete and experimentally allowed 
realization of the composite higgs theory. We will apply the ATLAS selec¬ 
tion to Monte Carlo generated samples of the heavy gluon signal, finding the 
region of the model parameter space, currently not excluded by experiments, 
which can reproduce the ATLAS excess. 

The paper is organized as follows: we present our benchmark model in 
sec. [ 2 j we discuss the experimental limits on the model in sec. |3j we apply 
the ATLAS selection to the composite Higgs signal in sec. [4] and present the 
results in sec. [5j We finally draw our conclusions and discuss the implication 
of our findings, especially for future experimental tests, in sec. [6j 


2 The model 


We evaluate the heavy gluon contribution to the signal examined by ATLAS 
in a concrete model. We consider the effective description of heavy gluon 
and vector-like quarks (vlq) interactions derived in 11,12 . This effective 


model can reproduce the low energy limit of compelling theories to solve the 
hierarchy problem and explain the electroweak symmetry breaking (EWSB), 


such as composite Higgs 13 or RS in a dual five dimensional picture |14|. 


The effective theory consists of two sectors, the elementary sector with ele¬ 
mentary particles and (analogous to) Standard Model (SM) gauge symme¬ 
tries and the composite sector, with composite particles, resulting from a 
new strong dynamics, which include the Higgs and new fermionic and vecto¬ 
rial resonances. The two sectors mix with each other. In particular, a new 
composite gluon associated with a SU(3) c ° mp gauge symmetry in the strong 
sector mixes with the elementary gluon of SU(3)f e . After diagonalizing the 
mixing, the physical states consist of the massless SM gluon and of a new 
heavy gluon G*. The G* interactions with SM particles and new vlqs from 
the composite sector are controlled by the parameter which describes the 
rotation of the elementary-composite mixing to the physical basis: 


tan 6*3 


n ele 
9 3 

gComp 


< 1 9s 


gf cos e 3 


comp • /) 

: g 3 p Sill 0 3 , 


( 1 ) 
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where g s is the QCD coupling and gf e ( g c 3 omp ) is the SU(3)f e (SU{3) c ° mp ) 
coupling. It is assumed that gl° mp > gf e . 


New vlqs also emerge from the composite sector. Based on the minimal 
SO(5)/SO(4:) coset for composite Higgs models 13 , which includes a cus¬ 


todial symmetry protection against large correction to the p parameter and 
to the ZbJ)L coupling (lO|, we consider composite fermions in fundamental 
representations of SO( 5). 5x representations decompose as (2, 2)x © (1, l)x 
under -S'0(4) x U{l)x ~ SU(2)l x SU{2) r x U(1)x ~ U(l)x is introduced to 
correctly reproduce the hypercharge, X = Y — T 3R . The composite fermion 
content we consider is thus 


Q = ( o T /3 ) = ( 3 ’ 2 ’ 2 ) 2 /3 T = (3,1, 1 ) 2/3 

\ B y (2) 

Q,= ( sli/3 B' ) = (3j 2j 2) - 1/3 ^ = (3 ’ V- 1 ' 3 ’ 

where we have specified the quantum numbers under SU(3) C x SU(2)l x 
SU(2) r x U(l)x- Similarly to the mixing for the vector fields, composite 
fermions also mix, via linear mass mixing terms |15), with the top and bottom 
quarks of the elementary sector. In particular, the doublet q\ = (tf e , bf e ) 
mixes with the (' T,B) SU{2)l doublet in the 52/3 and has a weak mix¬ 
ing with the (T', B') doublet in the 5 _i/3; t e ^ e , b \mix respectively with 
the electroweak singlets T and B. The Lagrangian of the mixing in the 
elementary-composite basis reads: 

C mix = -A L1 g£(T, B) - A L2 g| (T, B') - A R1 i R T - A R2 b R B + H.c. , (3) 


where A L2 X Li ~ A m ~ A R2 . The mixing leads to a scenario of partial 
compositeness of the top and bottom SM particles, which become superpo¬ 
sitions of their composite and elementary modes and acquire their masses 
through the interactions of their composite modes with the composite Higgs: 


nn t ~ Y,s l s r 


v 

71 


m b ~ Y*s 2 s bR 


V2 ’ 


(4) 


where Y* is the Yukawa coupling among composites, of 0( 1), and sl, s r , s b R 
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represent respectively the q R , the t R and the b R degree of compositeness [^} 
s 2 1 is a parameter associated with the q\ mixing with the (X 7 , B') SU ( 2 ) 
doublet in the 5 _i /3 second tower of resonance, which permits the genera¬ 
tion of the bottom mass. The weakness of this mixing justifies the smallness 
of the rrib/mt ratio and preserves the custodial symmetry protection to the 
Zb L b L coupling. We refer the reader to 111, 


12 


for more details on the model. 


In our study, for a reason which will appear clear in the next paragraphs, 
we are particularly interested in the bottom partners of the 5 _i /3 and we will 
consider the following simplifying assumptions on the vlq spectrum: 


Mg Mqi. The B decay modes depend on the details of the elec- 
troweak mixings among bottom fermions 11 16 . We thus decide to 


consider a decoupled electroweak singlet bottom-prime to reduce the 
model dependence of our analysis. 


• Mq Mqi. This is a simplifying conservative choice, since some of the 
heavy quarks in the Q, in particular the T and the X 2 / 3 , which decay 
~50% into Zt, can also contribute to the excess of events measured by 
ATLAS. 


As an effect of the custodial symmetry, which protects the Zbb coupling 
and demands an SU(2)l x SU(2) r bidoublet representation for the Q! com¬ 
posite fermions, a degenerate doublet of bottom partners, the B_ 1/3 and the 
B', with the same mass and Yukawa coupling (before the EWSB), appears 
in the spectrum. We can rotate the degenerate states to the new fields 16 


Bh - ^(-B-1/3 + B') 


B z =-(B„ 1/3 -B') 


(5) 


which will coincide with the mass eigenstates (up to a negligible correction 
for the B r , coming from electroweak mixing effects). The B R and the B z 

2 sl, sr, SbR are related to the mixing parameters in § by: 

Ail Ari Ar2 

Sl = — , = SR = = SbR — , = , 

\J ^ Li + Ml y / ^ 2 ri + Mj, \J ^- 2 r2 + 

where Mq , Mf, Mg are the Q, T and B masses before the elementary-composite mixing. 
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vlqs completely decay respectively into bh and into bZ (more details can be 

( 6 ) 


found in 16,17 and in the appendix): 

BR(B H -> hb ) = BR(B Z Zb) = 1. 


Notice that If we had a single SU(2) L doublet (T',B'), we would have had 
a single bottom-prime quark, the B', decaying 50% to Zb and 50% to hb, as 
explained in the appendix. In this case, the G* decays to B' pairs would have 
produced a mixed Zbhb final state 50% of the time. The custodial symmetry, 
which leads to the B H ( B z ) mass states fully decaying into hb (Zb) prevents 
the production of a Zbhb final state via the G* decays into bottom-prime 
paii'^0 This implies that the cross sections for the ZZbb and the hhbb final 
states are enhanced by a factor of two compared to the non-custodial cas^J 
The rest of the non-decoupled vlqs have the following decay patterns: 

BR(B _ 4/3 -»■ W~b) = BR(T' -t W + b) = 1 

BR(T ->• Wb) = 0.5 BR(T -)• Zt) = BR(T -> ht ) = 0.25 

In addition to the elementary-composite mixing, the EWSB induces a further 
mixing among same charge quarks. After the electroweak mixing diagonal- 
ization (the mass matrices can be found in jl8|), we finally arrive at the 
mass basis. We select the following set of mixing parameters, for which the 
bottom and top masses are correctly reproduced, and which give vlq (physi¬ 
cal) masses that fulfill the limits from direct searches for top-partners at the 
LHC: 


y* = 3 s L = 0.57 s R = 0.6 s bR = 0.3 (8) 

0The bare masses for T and for the Q! bidoublet have been fixed respectively 
at 880 GeV and 930 GeV. The physical masses are: 

M v = 945 GeV [912 GeV 


M b _ 4/3 = 930 GeV [912 GeV 


19 


19 


M f = 900 GeV [800 GeV 19,20 ] M Bh = 955 GeV [846 GeV 21 


M Bz = 930 GeV [700 GeV 22 ] 


3 except possible contributions from the G* decays into B pairs. 

4 assuming the same branching ratio for the heavy gluon decays to a generic pair of 
bottom-primes in the custodial and non-custodial cases. 

5 The remaining parameters, which we do not report in eq. ^ because are less relevant 
for the G* phenomenology, are: S 2 = S 3 = 0.03, S 4 = 0.05 (where S 3 , S 4 are small mixing 
parameters proportional to S 2 12 ). 
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We have indicated in parenthesis, for each vectror-likc quark, the strongest 
mass limit placed by LHC and a reference to the corresponding search. 


After having discussed the fermionic spectrum we now briefly examine the 
G* phenomenology. As anticipated, the heavy gluon interactions are ruled 
by the tan #3 parameter in ( 0 - The heavy gluon is essentially a composite 
particle, which thus interacts strongly with the composite modes. In partic¬ 
ular, the coupling of the interactions with elementary modes, thus with light 
quarks, is gs tan 83, that with composite modes, as the B z and Bh bottom 


partners, is gs cot 83. The G* branching ratios read 11 


BR(G* ->• jj) = -as tan 2 9 3 Mq* 

BR{G* —> tstL + bibs) = — r (sj cot 83 — cj tan# 3 ) Mq* 

o 

BR(G* t R t R ) = ( s 2 r cot 6» 3 - c 2 r tan 0 3 ) 2 Mq* 

a s 


BR(G* -> b R b R ) = — ( s 2 bR cot 0 3 - c 2 bR tan 0 3 ) M G * 


a s 


BR(G* -> Tt + Tt ) = —M g *~ 


s 2 r 2 
b R c R 


6 sin 2 83 cos 2 83 


M\ 

M 2 g * 


a s 


1 - - 


1 M% 1 M\ 


2 Ml, 


2 Ml, 


BR(G* TT) = ^ M g * 


cj j cot 8 3 — s R tan 83) + cot 2 83 

+6 (4 cot 2 9 ,-4 )Ti 


Ml 

Ml, 



Ml 
4—— 

Ml* 


inner -> Q'Q') - a fM G . cot 2 $ 3 (l + 2 dl^4 


K 

AP C , 


(9) 


where, c x = y/l — s 2 , x = L, R, bR and Q' denotes a generic resonance in 
the Q! bi-doublet: Q' = Bz , Bh, h>-i/4, T'. 

As an effect of the strong interaction with the heavy vlqs, G* decays predom¬ 
inantly into pairs of vlqs above the threshold 2 M v i q . In particular, in the G* 
mass region relevant to our analysis, we find 

BR(G* -> B Z B Z ) ~ 0.18 . (10) 
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Figure 1: G* production cross section at the yfs = 8 TeV LHC, for a value 
tan0 3 = 0.5. The cross section scales with tand 3 as tan 2 0 3 . 

The G* is produced by quark-antiquark annihilation (the gluon fusion pro¬ 
duction is suppressed by gauge invariance). The production cross section, 
due to the gs tand 3 coupling of the G* interactions with light quarks, de¬ 
pends quadratically on tan0 3 . We show in fig. [l] the G* production cross 
section at the y/s = 8 TeV LHC, for a value tan0 3 = 0.5. 

Summarizing the main results of this section, we have found that mini¬ 
mal composite Higgs theories with custodial symmetry protection to the p 
parameter and to the Zbb coupling and which can generate both the top and 
the bottom masses, predict a significant contribution to hhbb and ZZbb final 
states (enhanced by a factor of 2 compared to the non-custodial case) com¬ 
ing from heavy gluon decays to bottom partners. In particular, the process 
pp —> G* —> BzBz —» ZZbb will be the focus of our analysis. If one of the Z 
of the final state decays to letpons and the other to neutrinos, as shown in 
fig- n we can indeed have a significant contribution to the excess of events 
measured by ATLA^J 

6 The contribution from G* —> TT —>■ ZZtt to the ATLAS signal is found to be small, 
below the 5% of the total contribution, and will be thus neglected in our analysis. 
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Figure 2: Feynman diagram for the heavy gluon decays to Bz vector-like quark 
pairs, leading to a signature with two same flavor opposite charge leptons plus jets 
and missing energy. 


3 Current limits from direct searches for heavy 
gluon resonances 

In this section we discuss the portion of the G* mass vs coupling parame¬ 
ter space which has been excluded by current searches at the LHC. There 
have been no dedicated searches so far for heavy gluon decays to vector-like 
quarks and the only test at colliders for G* theories can be obtained from the 
searches for dijet and tt resonances. We will see that, since the G* branching 
fraction for decays to dijet or tt are small when the heavy gluon is above 
the threshold for decaying into pairs of vector-like quarks, the searches in 
the dijet and tt channels leave a large part of the model parameter space 
untested. 

The strongest constraints on tt resonances come from a recent CMS combined 
analysis of the searches in different final states 
in 


23 . We use the upper limits 


23 on the production cross section times tt branching ratio of a Kaluza- 


Klcin gluon to extract the exclusion region for the G* in the (Mg*, tan 6 * 3 ) 
plane. We find that searches in the tt channel can only exclude a portion of 
the parameter space at lower G* masses, below ~1.87 TeV. The small region 
of the parameter space excluded by tt searches which is of interest to our 
study is colored brown in our final plot of fig. [4j 

Searches for new particles in the dijet mass spectrum generally apply to 
narrow resonances with an intrinsic width much smaller than the experimen- 
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tal dijet mass resolution and cannot place limits on wide resonances like a 
heavy gluon decaying to top partners. CMS has however recently extended 
his search in the dijet channel to include an analysis for wide dijet reso- 
We use the results of this study to extract the limits on the 
In particular we use the upper limit, quoted in 24 


24 


nances 

G* parameter space 
on the production cross section times dijet branching ratio times acceptance 
(cr x BR x A) of a generic qq resonance with a width over mass ratio of 0.10. 
This is a conservative choice since, as we also show in fig. [4| the G* width 
over mass ratio is above 0.15 in most of the G* parameter space relevant for 
our study. Q As expected from the tan 2 0 3 scaling behavior of the G* pro¬ 
duction cross section, we find that the searches for dijet resonances exclude 
only a part of the G* parameter space at large (> 0.5) tan 6*3 values. The 
region of the G* parameter space excluded by dijet searches is colored grey 
in our final plot of fig. [4j 


4 Applying the ATLAS selection 

In this section we apply the main requirements of the ATLAS selection in |IJ 
to Monte Carlo generated samples of the process pp —> G* —> B Z B Z —> ZZbb 
in fig. [2] We generate the events for the heavy gluon signal at leading or¬ 
der with MADGRAPH 5 26 , after having implemented the model in sec. 
[2] by using feynrules 27|. We use the cetq611 PDF set |28| with renormal¬ 
ization and factorization scales fixed at the heavy gluon mass. The events 
are then passed to PYTHIA 6.4 |29j (with the default tune) for showering 
and hadronization. Jets are reconstructed with FAST JET 130 by an anti-kt 


algorithm with cone size R =0.4. In order to mimic detector effects we also 


More in details, we calculate the G* excluded regions by rescaling the a x BR x A 


values for Colorons/Axigluons ( C ) 25j shown in 24 according to 


0.86 x 


cr(G*) BR(G* -X qq) 
a(C) (C -> qq) 


= 0.86 x tan 2 83 x 


BR(G* -X qq) 
2/3 


the universal Coloron model considered by CMS 25 is indeed reproduced by the G 


model in the limit: tan0 3 = 1, Sr = Sr = s^r = 0 and decoupled top partners. The factor 
0.86 is included to take into account the limitations of the narrow width approximation, as 
suggested in 24 , and is taken from table 7 in [24]. The resulting values are then compared 


to the upper limits for qq resonances with T/M = 0.1. 
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apply a Gaussian smearing to the jet energy with: 


yg) ^, n , s 

E E + /E 


(11) 


where E is in GeV and C = 0.025, N = 1.7, S' 
turn is then rescaled by a factor E smeared / E. 


0.58 31 . The jet momen- 


Retracing the ATLAS analysis, we consider the following final state: two 
same flavor (electrons or muons) leptons with opposite charge plus at least 
two jets and missing energy: 

e + e~// j, + h~ + rij e t jets + E™ lss , rij et > 2 . (12) 

As a first step of the selection we apply the following isolation criteria and 
Pt requirements on the two final leading leptons (we apply the same cuts for 
electrons and muons): 


lh(^i, 2 )| < 2.4 Ai?(/i i 2 ,jet) > 0.3 AR(li,l 2 ) > 0.3 

p T h > 25 GeV p T k > 14 GeV . ^ ’ 

The jet separation requirement is applied to any “baseline” jet with px > 20 
GeV and \rj\ < 5. li (Z 2 ) denotes the leading (sub-leading) lepton. 

Signal jets must fulfill the conditions: 

|r/(jet)| < 2.5 pt jet > 35 GeV . (14) 

After these “acceptance” cuts, we get, at y/s = 8 TeV LHC with 20.3 
fb _1 , the number of G* signal events shown on the second column of table [lj 
The expected number of events is shown for several heavy gluon masses and 
for a fixed coupling tan # 3 = 0.5. 

The on-Z ATLAS search is focused on generalized gauge mediated supersym¬ 
metry breaking signatures, where pair-produced gluinos decay via neutralinos 
to a gravitino, which is the lightest supersymmetry particle, plus a Z boson. 
These GGM topologies are characterized by energetic final states with large 
missing energy associated with the gravitino. The ATLAS analysis thus 
applies the following set of cuts, aiming at selecting events with Z-peaked 
dilepton mass, large missing energy and energetic final states: 
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81 GeV < mu < 101 GeV E^ lss > 225 GeV 
H t > 600 GeV A0(jet 12 , E™ lss ) > 0.4 , 


(15) 


where H T is defined as: H T = prih) +PT(h) + ^2^^ PT(jet l ) , mu is the dilep¬ 
ton mass, E1^ , ' LSS denotes the transverse missing energy and A0(jet x 2 , E™ lss ) 
is the azimuthal opening angle between the missing energy and the leading 
or sub-leading jet; the restriction on the azimuthal jet -E™ lss separation is 
applied to reject events with jet mismeasurements contributing to large fake 
missing energy. 

We find that the ATLAS search strategy, with the main set of cuts in 
( [T5| , does not only apply to GGM topologies but also selects heavy gluon 
G* —s- BzBz —> ZZbb signals (fig. [ 2 ]), which are also characterized by a 
leptonically decaying Z , an energetic final state and large missing energy, 
produced by the decay to neutrinos of one of the two Z bosons in the final 
state. The H T and E™ zss distributions for the G* —> B Z B Z events in the 
channel (12) after the acceptance cuts in (13), (14) are shown in fig. [3] for 
several G* masses 0 

Applying the ATLAS selection (15) to the G* signal, we find the expected 
number of events, at the 8 TeV LHC with 20.3 fb 1 , indicated in tab. [TJ 
Columns 3-6 show the number of events passing the different steps of the 
ATLAS selection for several G* masses at a fixed coupling tan 0- A = 0.5. 


5 Results 


After the complete selection in (15) with 20.3 fb 1 of collected integrated 


luminosity, ATLAS observes an excess of events above the expected Standard 
Model background with a statistical significance of 3 standard deviations. 
The statistical significance is of 3 sigma in the electron channel and of 1.7 
sigma in the muon channel. The excess of events above background in the 
ee + /i/i channel, which can be read from tab. 7 in 0 , is: 


n 


= 18.4 ±3.2 


(16) 


8 We calculate E™ lss by Monte Carlo truth: we sum vectorially the transverse momen¬ 
tum of the neutrinos and those of soft (px < 20 GeV) or lost (|?/| >5) jets. 
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Figure 3: Normalized Ht (left plot) and E™ lss (right plot) distributions of G* —> 
BzBz events in the same flavor opposite charge dilepton plus jets plus missing 
energy channel, eq. (12), after the acceptance cuts of eq. (13), (14). Distributions 
are shown for G* masses of 1.9, 2.0, 2.1 TeV and tan $3 = 0.5. 


By applying the analysis explained in the previous section to G* —> BzBz 
events for different heavy gluon masses and couplings we can derive the 
(Mq* , tan 6(3) values which are able to reproduce the ATLAS results. The 
last column in tab. [l] indicates, for the different G* mass points, the tan 6 3 
values which give the excess of events in (16) measured by ATLAS. 


Our results are finally shown in fig. [4j The green band shows the ( Mq *, tan 63) 
values giving the excess of events in ( [16] ) within ±lcr from the central value, 
indicated by the black dashed curve. The grey upper region is excluded by 
searches for dijet resonances, calculated, as explained in sec. [3j from the 
CMS search |24|. The lower brown region is excluded by searches for tt res¬ 
onances and is derived (sec. [3]) from the CMS analysis in [23]. We see that 
heavy gluon decays to vector-like bottom quarks can explain the ATLAS 
finding in a parameter space region so far untested and thus not-excluded by 
LHC searches. We find that the ATLAS excess of events can be produced 
by G* —y BzBz events for G* masses roughly in a range 1.87 - 2.15 TeV and 
for gs tan 63 couplings within (0.3 - 0.65)<7 s. 


Finally, we show in Lg.s 5|6|7 the kinematic distributions of G* —>■ BzB- z 


events, to be compared with the distributions of the observed events shown by 
ATLAS in H] - which we report for the reader convenience on the right side 
of hg.s [5|6|7[ Since the highest statistical significance is found in the electron 
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M G * (GeV) 
[tan #3 = 0.5] 

Tljet 

> 2 

mi 

[81 , 101 ] 

TT'miss 

JUJrj-i 

>225 

Ht 

>600 

A^jet 1 j2 ,E™ ss ) 
>0.4 

tan 6*3 

[n exc = 18.4 ± 3.2] 

1870 

68 

60 

52 

51 

46 

0.29 ± 0.03 

1900 

56 

49 

42 

42 

38 

0.33 ± 0.03 

1950 

38 

33 

28 

28 

25 

0.41 ± 0.04 

2000 

30 

26 

23 

23 

20 

0.47 ± 0.05 

2100 

21 

19 

16 

16 

14 

0.58 ± 0.07 

2150 

18 

16 

13 

13 

12 

0.65 ± 0.09 

2200 

16 

14 

12 

12 

11 

0.74 ± 0.14 


Tabic 1: Number of events at the 8 TeV LHC with 20.3 fb _1 for the G* —> BzBz 
signal in the same flavor opposite charge dilepton channel which pass the different 
steps of the ATLAS selection JTj (mu, E™ lss and Ht cuts are in GeV). We show 
the results for several G* masses at a fixed coupling tan $3 = 0.5. The last column 
indicates, for each G* mass point, the tan $3 value which gives the excess of events 
measured by ATLAS |T]. 


channel, we just show the distribution in the e + e~ channel. We find that 
G* —> BzBz event distributions are compatible with the jet multiplicity, the 
H t , the E™' lss , the mu and the A(/>(jet 12 , E™ lss ) distributions of the ATLAS 
observed events □with few exceptions: the number of observed events on the 
first bin in the Ht distribution and on the first bin in the E™ lss distribution 
are slightly above the expected number of G* — * BzBz plus background 
events^] In general, we find that the Ht, the E™ lss and the ?n« distributions 
for G* events are similar to those for the GGM model with tan f3 = 1.5, gluino 
mass of 900 GeV and p = 600 GeV, shown by ATLAS [I]. Most notably, 
the jet multiplicity distribution of G* —> BzBz events looks closer to that 
of the observed events compared to the jet multiplicity distributions for the 
two examples of GGM models considered by ATLAS, which tend to predict 
a number of signal jets larger than the one observed. 

9 Obviously, G* signal events have to be summed to the expected background events in 
order to be compared with data. 

10 Here we just comment qualitatively on our findings for the event distributions. A 
detailed analysis of the statistical compatibility between the G* —> BzBz and the observed 
distribution of events is beyond the scope of this paper and we leave it to future studies. 
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Figure 4: (Mg* , tan 6 h) values reproducing the ATLAS |1| excess of events. The 
green band shows the values giving the excess of events within Tier from the central 
value (black dashed curve). The grey upper region is excluded by searches for dijet 
resonances - it is derived from the CMS limits on qq resonances with T/M = 0.1 
; the brown region on the left is excluded by searches for tt resonances |23j. 
Dotted curves show (Mg*, tan$ 3 ) regions where T(G*)/Mg*= 0 . 10 , 0.15, 0 . 20 . 

6 Conclusions and discussions 

In this study we have found that the 3 sigma excess recently measured by 
ATLAS in events with same flavor opposite sign lepton pairs, peaked on the 
Z mass, plus jets and large missing energy can be interpreted in composite 
Higgs models as a result of heavy gluon decays to vector-like quarks. We 
have analyzed a concrete model where, as an effect of a custodial symmetry 
protection to the Zbb coupling, bottom partner vector-like quarks, produced 
in pairs from the heavy gluon decays, contribute significantly to hhbb and 
ZZbb final states. This latter ZZbb signature, in particular, if one of the Z 
decays leptonically and the other to neutrinos, can explain the ATLAS ex¬ 
cess. Our results are summarized by fig. [4] where we show that the ATLAS 
excess can be reproduced by the composite Higgs model, in an experimen¬ 
tally allowed region of the parameter space, for heavy gluon masses roughly 
in a range 1.87 - 2.15 TeV and for heavy gluon couplings to light quarks 
within ~ (0.3 — 0.65)g,s. 

An immediate reaction to this finding is looking for confirmations of the com¬ 
posite Higgs interpretation of the ATLAS results, especially to distinguish 
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Figure 5: Kinematic distributions of events (only for the ee channel): jet multi¬ 
plicity in the upper plots and Ht in the lower plots. The distributions are obtained 
after all of the cuts in (15) have been applied. The Plots on the left refer to the 
pp —>• G* —» BzBz —> ZZbb signal for different heavy gluon masses (1.9 TeV, 
2.0 TeV, 2.1 TeV) and for the tan @3 values reproducing (the central value of) the 
excess of events measured by ATLAS l] in the ee + pp channel. The Plots on the 
right are taken from ATLAS 111 and show the jet multiplicity and Ht distributions 
of the observed events, of the estimated backgrounds and of two examples of GGM 
models with tan B= 1.5. 
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Figure 6: Kinematic distributions of events (only for the ee channel): dilepton 
mass in the upper plots and missing transverse energy in the lower plots. The 
distributions are obtained after all of the cuts in (15) have been applied. The 
Plots on the left refer to the pp —> G* —» BzBz —> ZZbb signal for different heavy 
gluon masses (1.9 TeV, 2.0 TeV, 2.1 TeV) and for the tan 63 values reproducing 
(the central value of) the excess of events measured by ATLAS 1 in the ee + pp, 
channel. The Plots on the right are taken from ATLAS (T| and show the mu and 
EJp iss distributions of the observed events, of the estimated backgrounds and of 
two examples of GGM models with tan/3= 1.5. 
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Figure 7: Kinematic distributions of events (only for the ee channel): azimuthal 
opening angle between the missing energy and the leading (upper plots) or sub¬ 
leading (lower plots) jet. The distributions are obtained after all of the cuts in 
(15), except the requirement A<h(jet 12 , E™ lss ) >0.4, have been applied. The Plots 
on the left refer to the pp —> G* —> BzBz —>• ZZbb signal for different heavy gluon 
masses (1.9 TeV, 2.0 TeV, 2.1 TeV) and for the tan $3 values reproducing (the 
central value of) the excess of events measured by ATLAS jTj in the ee+pp channel. 
The Plots on the right are taken from ATLAS |[lj and show the A4>(jet 1;2 , E™ lss ) 
distributions of the observed events, of the estimated backgrounds and of two 
examples of GGM models with tan /?= 1.5 
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Figure 8 : Event distributions (with 20.3 fb ^ 1 at the 8 TeV LHC) for the G* —> 
BzBz signal with parameters reproducing the ATLAS excess. The vlq mass is 
fixed at Mb z = 930 GeV. Left plot: invariant mass of the two leptons plus the 
leading jet, for Mg* = 1.9 TeV (tan 63 = 0.33). Right plot: invariant mass of 
the two leptons plus all of the signal jets, for Mq* = 1.9, 2 . 0 , 2.1 TeV (tan #3 = 
0.33,0.47,0.58). 

the composite Higgs hypothesis from alternative supersymmetrical explana¬ 
tions. The s-channel exchange of the heavy gluon and the two vector-like 
quark resonances, one of which can be fully reconstructed, are the distinctive 
features of the composite Higgs signal. Indeed, a simple test of the composite 
Higgs interpretation could be easily realized by analyzing, for example, the 
event distributions of the following observables: the invariant mass of the two 
leptons plus the leading jet, which shows a peak around the vector-like quark 
mass, as shown on the left plot of fig. [ 8 j or the invariant mass of all of the 
observed objects, leptons and signal jets, of the final state, which presents a 
kinematic edgfp] at high mass, as a result of the heavy gluon exchange; the 
corresponding plot is shown in fig. [ 8 | 

If the ATLAS excess is really due to a composite Higgs/RS theory, further 
evidences should manifest soon at the upcoming LHC run. The bottom and 
top partners of the model should be indeed observed in searches for pair 
production of vector-like quarks. An evidence for the electroweak singlet T 

11 because of the missing energy, it is not possible to fully reconstruct the heavy gluon 
resonance. 
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could also appear in searches for top-partner clectroweak single production; 
dedicated analyses have been performed in |32--j"34p In general, the results 
of this study reinforce the importance of focusing on searches for vector res¬ 
onances which include the decays to vector-like quarks 11,35-37 . A large 


portion of the parameter space for composite Higgs models is indeed not 
accessible to standard search channels, as those in the dijet or tt final state. 
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Appendix 


In this appendix we demonstrate that BR(Bz — > Zb) = BR(Br —> hb) = 1 . 
This result can be easily obtained from the Yukawa Lagrangian, which in the 
elementary-composite basis reads |l2|: 

C YUK — Y*Tr[Q'H]B -I-, (17) 


where we have written only the terms relevant for our purpose. R is the 
Higgs matrix that, written in terms of the Higgs and Goldstone bosons, is 


n 


h — iz a/2w + 
— y/2 ui~ h + iz 


( 1 . 2 , 2 ),/,. 


(18) 


After the diagonalization of the elementary-composite mixing, which leads 
in particular to the new B eigenstate: 


Br — c bR B c ° mp + s bR,bR e 


SbR 

CbR 


A 


R2 


M 


B 


bR — —SbRB R mp + CbRbft C 

PI the Yukawa Lagrangian becomes 

' A = -Us SR [-§-1/3,, (h - iz) 6r + B’ l (h + iz) b R ] + 


(19) 


( 20 ) 


12 with the choice of parameters in eq. (|8j) , the electroweak coupling of the T is Aj. ~ 1.4. 
For such a coupling the T should be discovered with 100 fb -1 at the 14 TeV LHC for masses 


up to ~ 1 TeV 32 . 

13 here Mr denotes the bare B mass, before the elementary-composite and the elec¬ 
troweak mixings. 


20 













After the field rotation in (5), B H ( Z ) = 
C YUK = -y/2Y,s m [B Hl 


^[B- 1/3 + (-)B'], we finally have 
hbn — iBz L zbft\ + • • • . ( 21 ) 


According to the equivalence theorem, which can be safely applied in our 
case, since we work in a regime M Rh ~ v, the above equation 

implies 

BR(B Z Z L b) = BR(B h -y hb) = 1 . (22) 


We can also notice that, in a non-custodial scenario, with only the B’ bottom 
partner, eq. (20) simply reduces to: 


C YUK = -Y*s bR [B' L hb R + iB> L zb R \ + • • • (23) 

and thus the B' decay pattern is: BR(B r —» hb) = BR(B r —» Z R b) = 0.5. 
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